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We investigate the state space of bipartite qutrits. For states which are locally 
CN ■ maximally mixed we obtain an analog of the "magic" tetrahedron for bipartite 

qubits — a magic simplex W. This is obtained via the Weyl group which is a kind of 
"quantization" of classical phase space. We analyze how this simplex W is embed- 
\ ded in the whole state space of two qutrits and discuss symmetries and equivalences 

inside the simplex W. Because we are explicitly able to construct optimal entangle- 
ment witnesses we obtain the border between separable and entangled states. With 
our method we find also the total area of bound entangled states of the parameter 
subspace under intervestigation. Our considerations can also be applied to higher 
\£> ' dimensions. 
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d • I. INTRODUCTION 

In 1935 Erwin Schrodinger stated already that "entanglement is the quintessence of the 
quantum theory". The late discoveries and developments in many distinct branches of 
physics show its immense validity. It is the basis for quantum cryptography, quantum 
teleportation and maybe if realizable quantum computation. It has also triggered a new 
field: quantum information. 

The main problem for composite systems is how to find out if a given state is separable or 
not and thus to characterize the border between separability and entanglement. While we 
have for the simplest composite system — two two-level systems (2x2 systems) or bipartite 
qubits — a necessary and sufficient criterion for separability, the Peres criterion, it is for 
higher dimensions only a necessary criterion (except 2x3). The criterion states that every 
separable density matrix is mapped into a positive semidefinite matrix by partial transposi- 
tion (PT), i.e., by a transposition on one of the subsystems. The reason why it fails for higher 
dimensions is that for these systems a completely positive map cannot be characterized by 
transposition alone and moreover these systems show more aspects of entanglement. 
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It is obvious that the knowledge of the state space is the key ingredients to understand 
entanglement and therefore for developing and optimizing new applications. Moreover it 
will help in understanding the relation of different entanglement measures. 
In this paper we focus on bipartite qutrits (3x3 systems). We concentrate on the set of 
locally mixed states with a quasiclassical structure and construct a geometrical picture of 
the state space. The quasiclassical structure fits also exactly into the conditions needed for 
teleportation and dense coding, e.g. Refs. [3, 0, 0|- While these sets of states have been 
noted already in Ref. (HHI, only little is known about its structure concerning entanglement, 
witnesses, PPT (positive partial transposition) and possible bound entanglement. 
For two qubits four orthogonal Bell states can be used to decompose every locally mixed 
state and a geometric picture can be drawn. In such a geometrical approach the Hilbert- 
Schmidt metric defines a natural metric on the state space, e.g. Ref. (a III- Via diagonalizing 
every locally mixed state can then be described by three real parameters which can be used 
to identify the density matrix by a point in a 3-dimensional real space. The positivity 
condition forms a tetrahedron with the Bell states at the corners and the totally mixed 
state, the trace state, in the origin. Via reflection one obtains another tetrahedron with 
reflected Bell states at the corners, see e.g. Ref. The intersection of both simpleces 

gives an octahedron where all points inside and at the border represent separable states 
because they are the only ones invariant under reflection and thus positive under PT. While 
for qubits this characterizes the separable set of locally mixed states fully we show in this 
paper that the analogue to the octahedron for qutrits is not quite that simple and in addition 
not all locally maximally mixed states can be imbedded. 

The simplex for bipartite qutrits lives in a 9-dimensional Euclidean space where the borders 
are given by the positivity condition of density matrices. We construct two polytopes and 
prove that they are an inner (kernel polytope) and an outer fence (enclosure polytope) to 
the border of separability. The boundary achieved by taking the set of all states which are 
positive under PT has not only linear faces and corners but also curved parts. Then we 
explicitly show how to construct optimal witnesses and apply them to certain states and 
show that there are regions where there is bound entanglement, i.e. entanglement which 
cannot be distilled by local operation and classical communication (LOCC). 
The paper is organized as follows, we present first the construction of the set of states we 
are analyzing, the magic simplex W. Then we discuss how the set is embedded in the whole 
set of states. We proceed with analyzing the rich structure of symmetries inside W: The 
symmetry of a discrete classical phase space. Hereupon we focus on describing the boundary 
of separability by calculating optimal witnesses. The optimization is done analytically and 
also numerically. Further we added an appendix for more details. 

Many of our considerations can be extended to pairs of qudits. In order to be as concrete 
as possible we postpone generalizations to higher dimensions and more abstract analyzes to 
a following companion paper. 

II. THE CONSTRUCTION OF THE MAGIC SIMPLEX W 

Throughout this paper we focus on two parties with 3 degrees of freedom e.g. "qutrits". 
Take any maximally entangled pure state vector in the Hilbert space C 3 <8> C 3 for defining a 
"Bell type state". Denote this vector as fi ,o- Choose the bases {|0), |1), |2)} in each factor 
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such that 

fi o,o = 73^|s>®|s>- (!) 

s 

On the first factor in the tensorial product - the side of Alice - we consider actions of the 
Weyl operators. They are defined by 

W k/ \s) = w k ^\s-e), (2) 
w = e 2 ™ /3 . (3) 

Throughout this paper the letters {s, t,j, k, £, m, n,p, q} denote the numbers 0, 1, 2. Calcu- 
lations with them are to be understood as "modulo 3". So "1 + 2"= 0, "2 x 2"= 1, 
"-1"=2, etc. 

The transformations which we consider take place on the first factor, on the side of Alice. 
Bob's side is in our definitions inert. This is not really an asymmetry. For the Bell state f2o,o 
every action of an operator A on the side of Alice is equivalent to the action of a certain A on 
the side of Bob. Concerning the Weyl operators the equivalent action is W k j = w~ 2M Wk-i- 
Changing the roles of Alice and Bob i.e. the flip transformation is therefore equivalent to 
a local reflection combined with a phase factor, but with no change of the total set of the 
produced states. The phase factors will disappear in the projection operators to be defined 
in equation (jEJ), and the reflection is one of the symmetries studied in Sect. IIVI 
The actions of the Weyl operators - we simplify the notation and write W^lfi) meaning 
(Wk t t <S> - produce on the whole nine Bell type state vectors 

tt k/ = Wfc,^o,o • (4) 
The Weyl operators obey the Weyl relations 

Wj^Wk.m = w ki W j+k ,£ +m , (5) 

Kt = = wMw -k,-e , (6) 

W ,o = 1 • (7) 

We remark that the Weyl operators and the unitary group which they form appear sometimes 
in disguise, under names like "generalized spin operators", "Pauli group" and "Heisenberg 
group", Refs. 0,0. 

The original use of the Weyl operators, in the chapter "Quantenkinematik als Abelsche 
Drehungsgr upp e" of Ref. 1101 was the "quantization" of classical kinematics. (Both continu- 
ous and discrete groups have their appearance there.) In the appendix we present a physics 
model for the bipartite system of qutrits, which may help to visualize the ideas, concerning 
the interplay of quasiclassical and quantum structures. 

The set of index pairs (k, £) is the discrete classical phase space; I denotes the values 
for the coordinate in "x-space" , k the values of the "momentum" , see also Fig. ^ 
To each point in this space is associated a projection operator 

Pk,e = \^k,e) (Qk,e\ ■ (8) 

This projection operator is the density matrix for a Bell type state. The mixtures of these 
pure states form our object of interest, the magic simplex 

W = { J]c M P M | c M >0, ^c M = l }, (9) 
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with the nine pure states P^/ at the corners. As a geometrical object it is located in an 
8-dimensional hyperplane of the 9-dimensional Euclidean space {A = ^ dk,ePk,e I a k,e £ 
K.}, equipped with a distance relation A/Tr(A — B) 2 . Specifying the origin A = in this 
Euclidean space, it is also equipped with a norm, the Hilbert-Schmidt norm VA 2 , and the 
inner product Tr(AB) = ^2 ak,i bk,e- 

The geometric symmetry of this simplex for the qutrits is larger than the symmetry which 
is related to the underlying algebraic relations. The later one is equal to the symmetry of 
the classical phase space. This is studied in Sect. II V I 



III. HOW IS W EMBEDDED IN THE SET OF STATES? 

W contains only states which are locally maximally mixed, i.e. every partial trace gives 
the unit matrix times the normalizing constant. Further it contains the maximal possible 
number of mutually orthogonal pure states. While this characterization is sufficient for 
qubits, defining the magical tetrahedron or any locally unitary transform of it, this is not 
so for the qutrits. More explicitly for qubits every locally maximally mixed state can be 
embedded into a magical tetrahedron, while for the qutrits we observe: 

1. There exist locally maximally mixed states, which cannot be diagonalized with maxi- 
mally entangled pure states, the Bell type states. 

2. Even if such a maximally mixed state is decomposable into orthogonal Bell type states, 
it may be inequivalent to any of the states in W. 

3. There are maximal sets of nine mutually orthogonal Bell states, which do not build 
an equivalent of W. 

Examples are presented in the appendix. 

We remark that there are other ways to characterize the density matrices, by expanding 
them into products of operators which are a basis for the space of matrices. The use of 
products of Weyl operators in Ref. [5[ is closely related to the construction in this paper. 



And it is the analogue to the use of products of Pauli matrices, e.g. Ref. , considered as 
forming a group. Another method has been tried, considering the analogue of Pauli matrices 
as generators of SU(2). This leads to using the Gell-Mann matrices, e.g. Ref. generators 
of SU {3). 

There are several ways how to characterize a special unitary equivalent of one of the versions 
of W. One is already given by the construction: Choose one of the Bell states, and choose 
some basis on one side. Another way would be a choice of fixed special Bell states which 
have to be represented. There is a precise statement about the restrictions and the freedom 
to do this: 

THEOREM 1 Every pair of mutually orthogonal Bell states can be embedded into a version 
ofW. Such a pair fixes the appearance of a certain third Bell state. A fourth state can then 
be embedded, if it is orthogonal to the first three. Then, with four different Bell states, all 
elements of W are fixed. 
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PROOF Choose one vector out of the pair as f2 ,o- Take a Schmidt decompositions of this 
vector and of $ the second one with the same basis on Bob's side: 

s s 

Consider the unitary operator U = J2 S \i/j s ){(j) s \, acting on the first factor. Orthogonality of 
the Bell states implies 

3(Q ,o|$> =Tr[/ = 0. 

This is possible only if the eigenvalues of U are the three numbers {e lS w k } with some 
common phase factor 5. Now let {\s)} be the eigenvectors of U, and fix Wi t o = e~ U . This 
implies Pi j0 = |$)($|- Note, that there are still three phase factors not fixed, one for each 
basis vector. Nevertheless, the projector P 2 ,o = WipPipW\ is defined unambiguously. 
The fourth Bell state has a Schmidt decomposition ^gX^IXs) ® I s )- Observe that 
the orthogonality to the first three states implies ^2i S w ks (s\x s ) = for each k, and so 
( s \Xs) — for each s. Together with the orthogonality of the Xs, this implies that either 
\Xs) = e ir, ^\s + 1), or \xs) = e iV ^\s + 2). Now fix the phases for each \s), such that all 
rj(s) = rj, and the fourth Bell state is implemented as either P 0)1 or as P 0)2 . So all the 
ingredients for the construction of W are fixed. □ 

That the special choice of the positions in the phase space makes no difference for 
the total set of elements is made clear by consideration of the symmetries inside W. 

IV. SYMMETRIES AND EQUIVALENCES INSIDE W 

We consider linear symmetry operations mapping W onto itself that can be implemented 
by local transformations of the Hilbert space. So separability remains unchanged. We show 
that the transformations of W can be considered as transformations of the quasi classical 
discrete phase space. 

Letting the Weyl operators act on Alice's side gives the phase space translations: 

%n,n '■ Pk,l l_> Pfc+m,£+n = W m>n P k ^W m ^ n . (10) 

The action is a discrete "Galilei transformation" . The quantization, expressed in the phase 
factors of the Weyl relations, disappears due to the combined action of W m ^ n and its adjoint. 
In the appendix we present the Weyl operators in matrix form, and also their relations to 
the phase space. As usual, we take the "x-coordinate" t as horizontal and the "momentum 
coordinate" k as vertical. 

For all the other special operations Po,o stays fixed. This is no general restriction, since 
translations may shift each of the points in phase space to the origin. Now we have to use 
the help of operators acting on Bob's side. For every linear operator A on Alice's side there 
exists an operator A acting on the other party, such that 

A\n 0:0 ) = A\n 0:0 ) . (ii) 

They are related through transposition in the preferred basis: 

(s\A\t) = (t\A\s). 
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FIG. 1: Here we plotted the points Pi yk of the discrete classical phase space. I denotes the values 
of the x coordinate and runs from to 2 and k "quantizes" the momentum and runs also from 
to 2. From one fixed point, e.g. Po,o, all possible lines are drawn. Thus the phase space can be 
divided into four bundles where each bundle consists of three parallel lines. In Sect. IIVI we show 
that transformations inside the simplex W are equivalent to transformations in this phase space 
and that the lines are all equivalent in the sense that each line may be transformed into any other 
one. This enables us to study the geometry of separability in the magic simplex W. 

So, for every local unitary U 

Utf\tt Q ,o) = |fio,o> ■ (12) 

Thus every unitary transformation of the Weyl operators can be lifted to a unitary trans- 
formation of W: 

W m , n = e^UWkjrf, (13) 
=> P m , n = UU ] P k/ U ] U. (14) 

This follows in detail from 

utfp k4 rtu = uu^w k/ p 0fi wl e uu^ = uw k4 frp 0fi uwi/P 

= uWwrfutfPQflrtuuwhrt = w m>n p 0j0 wi !n . (is) 

First consider 

f 

It effects the quarter rotation of phase space (counter clock-wise) 

K: P k/ ^Pe,- k . (16) 

Then consider 

U v : 10)- |0), |2>-u> 2 |2>, 
it lifts to the vertical shear 

V: P k ,^P w . (17) 
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Combined action perfects the horizontal shear 

n = n- 1 v- 1 n-. P M ^P M+fc . (18) 

Now, for the following reflections we have to consider anti-unitary transformations of the 
Hilbert space - unless we want to use the flip, the exchange of Alice's and Bob's side. The 
simplest one in our preferred basis is the vertical reflection 

S: P M ^P_ M . (19) 

It is implemented by complex conjugations 

E a *i s >^E^i s > ( 2 °) 

s s 

on both factors. Their tensorial product gives global complex conjugation 

C: Ea s , t |s)®|*)^E< t |s)®|*). (21) 



s,t s,t 



Obviously 



C- 1 = C, CW M |tf) = W_ M C|tt), P-k,t = CP k/ C. 



We notice that this anti-unitary transformation is also compatible with the Weyl relations. 
Since it acts globally, on both factors, it is positivity, separability and PPT preserving. All 
the structural properties of W which are of interest are symmetric under vertical reflection. 
Hence they are also symmetric under combined action with other transformations, which 
give new kinds of reflections, as horizontal reflection 

nsn- 1 : p m ^ p ki _ t (22) 

and diagonal reflection 

US: P k/ ^P e , k . (23) 

All the transformations are "linear" or "affine" mappings of the phase space. Phase space 
lines ((k,£), (k + m,£ + m), (k + 2m,£ + 2m)) are mapped onto lines. Note that the sequence 
of the three points in a line can be rearranged in any way. With the right relabelling of the 
indices modulo three one gets again the special prescribed form. 
Let us collect our results and state the following theorem: 

THEOREM 2 The group of symmetry transformations ofW is equal to the group of affine 
transformations of the quasi classical phase space which is formed by the indices, 



where mq — pn ^ 0, with all calculations done with integers modulo 3. For mq — pn = 1 the 
transformation of the Hilbert space is unitary, for mq — pn = — 1 it is anti-unitary. 
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PROOF Pure phase space translation by (J,r) is the second part, combined with the 
unit matrix, m = q = I, n = p = 0. The generating elements for the first part of 
transformations and the corresponding matrices are (in this Proof we use "—1" for "2") 

Quarter rotation: 1Z <-> | J j ' 

Vertical shear: V <-» f o 1 y ' 

Vertical reflection: S «-> ( ^ ^ J . (27) 

All the invertible matrices can be generated. This can be seen first by looking at the 
numbers of zeroes. Maximally two zeroes, in relative diagonal positions, are possible, as in 
1Z and S. One zero is possible, as in V. There may be different positions of the zeroes, 
but they can be moved by the diagonal reflection, applied from the left and/or from the 
right. The case with no zero in the matrix is represented by VH. Finally there are different 
distributions of signs, but they cannot be changed individually, since this would not give 
invertible matrices. The signs can be changed pairwise, for each column or row, by the 
vertical reflection S, and by TZSTZ' 1 , applied from the left and/or from the right. □ 

The group structure of the combined transformations can be written in matrix nota- 
tion: 

k\ I m n j \ I k 

lU L g r f | . (28) 
lj \o l) \l 

The lines in the discrete phase space play also an important role in connection with the 
Mutually Unbiased Bases, see Refs. 0, [l2[ EH . In Fig. [T] we visualize all possible lines 
for one phase space point. Thus we have for the whole phase space four bundles — called 
"striations" or "pencils" in Refs. 05 E3 > respectively — , each one with three parallel lines. 

This makes 12 special sets out of „ =84 subsets with three points of the phase space. 



v 3 / 

The lines are all equivalent in the sense that each line may be transformed into any other 
one. More general, we have 

THEOREM 3 In the classes of subsets of phase space points, there is just one equivalence 
class of single points, one of pairs, two classes of triples and two of quadruples. The equiva- 
lence relations are moreover valid also for the complementary sets with five to eight points. 
Inside each pair and inside each triple, there is total symmetry under permutations. 

PROOF We move the subsets to special chosen places in phase space. Consider one point of 
the subset after the other, in any order. Translation brings the first one to the origin (0, 0). 
The shear transformations bring the second one to (1,0). If these two are part of a line, the 
third one has then been moved automatically due to the linearity of the transformations, 
together with the first two, to the point (2, 0), completing this vertical line. If the third (or 
fourth) point is not in a line with the first pair, it is movable with horizontal reflection and 
vertical shear to the place (0, 1), without changing the arrangements of the line with i = 0. 
In the case of four points, with the vertical line at £ = not yet completed, we have several 
cases: If the fourth point is either at (2, 2) or at (0, 2), it completes another line, and we can 
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start again, moving this line to the preferred vertical one, and the remaining point as done 
above. In case the fourth point is not yet at (1, 1), where we want to place it, completing 
a square, it is either at (2, 1) or at (1, 2), and it can be moved by shear, together with one 
of the others, to form the preferred square. These are the cases, where no complete line is 
contained in the subset of four. 

The inner symmetries of pairs and triples are now implicitly proven, since the sequence of 
moving their points can be chosen arbitrarily. □ 



V. THE GEOMETRY OF SEPARABILITY 

We now proceed to the question of separability. We start with a rough analyzes of an inner 
and outer fence in W. Then we concentrate on the border given by PPT. In particular we 
show that the test for positivity under PT reduces to a check for positivity of a 3 x 3 matrix. 
As an example we study the entanglement of mixtures of the total mixed state and two 
orthogonal Bell type states. We then explicitly describe the construction of witnesses and 
analyze the strategy to optimize them. We apply our method to the above state and find 
for some mixtures bound entanglement. As an further example we discuss a density matrix 
which is a mixture of the total mixed state and three orthogonal Bell states where two of 
them are equality weighted. 



A. Two polytopes as inner and outer fences for separability 

The most mixed separable state, with density matrix u> = ~1, lies at the center of W. 

W= g-E P M, (29) 

k,£ 

since the Qk,e form an orthonormal basis. The separable states with the largest distance to 
the center are defined by the lines in the phase space. 

THEOREM 4 The twelve outermost separable states in W have the density matrices 

Pline = g Pk,e ■ (30) 

(k,£) e line 

PROOF This is a special case of the more abstract general statement in Ref. jl| equation 
(36). For a more concrete demonstration, consider the vertical phase space line {(k, £)} = 
{(0,0), (1,0), (2,0)}: 



lj2w k ^\s,s){t,t\, (31) 

s,t 

where we now write \s, t) for \s) ® \t). With J2k wk ^ s ~^ = 38 s ,t one gets 



Pk,e=o g 



k 

obviously a separable state. 



Pline = ^^2 P k,0 = ^^2\S,S)(S,S\, (32) 
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It lies in each one of the three hyperplanes B p0 in our Euclidean space, defined by 

B P ,q = {(ckd\cp*= 3}- ( 33 ) 

B Pt o intersects the line of isotropic states (1 — a)u + aP Pi o, exactly at the border between 
separable and entangled states at a = \, see also Ref. jj]. This hyperplane is therefore the 
proper entanglement witness, reduced to our subspace of hermitean matrices. Each state 
outside is entangled, and it is only the center of the triangle with the P Pi o at the corners 
which is a separable state. 

By the equivalence relations stated in Theorem El all these considerations are valid for 
each one of the twelve phase space lines. Now the witness hyperplanes intersect also at the 
centers of the other 72 (= 84 — 12) triangular faces, but the states there are not separable. 
This is easily checked by showing that they are not PPT. This is done explicitly in the next 
Sect. EE □ 

The nine pairs of hyperplanes 

Bp,q — {{ c k,i) I c p,q 

and A Ihq = {(c k/ ) \ c p , q 
enclose all the separable states of W. They define the 

enclosure polytope {(c/c/) | all c VA G [0,-]}. (36) 

3 

It has the same geometric symmetry as the simplex W, which has nine corners. Intersections 
of the -Bp )(J -hyperplanes in triples give 84 vertices. In each of the B PA there lie 28 of these 
vertices, the other 56 lie in A Ptq . 

Of these 84 vertices, only twelve are separable states. All the convex combinations of these 
twelve are again separable states. They form the 

kernel polytope {p = J~] Kpiinea I K > 0, J^A a = l}. (37) 

lines a 

It has twelve vertices, which are the pu ne - Each hyperplane B pq contains four of them. 
The other eight are in A p jq and define a full seven-dimensional convex body. This may be 
compared to the four triangular faces of the qubit-octahedron lying inside the triangles of the 
magical tetrahedron. The other four, out of all eight, lie in witness planes, see Refs. IllL Il4l . 
Here is one more of the many differences between qubits and qutrits (compare with Ref. jl5|): 
We do not have the geometric rotation-reflection symmetry between the bordering planes. 
In the witness-hyperplane B pq there is only a three-dimensional face (a tetrahedron) with 
four vertices of the kernel polytope. 

That the bordering face in A of the kernel polytope is seven-dimensional can be seen 
by looking in the Euclidean space of hermitean matrices at the eight vectors v p ^ q , with 
(p, q) 7^ (0,0). Such a vector v p , q is defined as pointing from |(1 — Po,o), the center of the 
face of W in v4o,o, to P Pj(? , one of the eight vertices of the seven- dimensional simplicial face of 
W. The center of this face of W is also the center of the eight vertices in A§ § of the kernel 
polytope. The vector is now representable as a linear combination of the pu ne , phase space 



1 



3 
0} 



} 



(34) 
(35) 
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lines through (p,q) not containing (0,0). By equivalences and symmetries it is sufficient to 
demonstrate this for one example: 

2^2,2 = -[(Plfl + P 2 ,l + P0,2) + {Po,l + Pl,2 + P2,0)\ + (38) 
-~[(P ll0 + Pi,! + Pl, 2 ) + (P 2 ,0 + Pl,l + P , 2 ) + (P ,l + Pl,l + P 2 ,l)] + 

+ ^[(P 2 ,0 + P 2 ,l + P 2)2 ) + (Pl,0 + Po,l + P 2)2 ) + (P ,2 + Pl,2 + P 2)2 )] • 

This consideration is, by equivalence, valid for the maximal face in any A Ptq . 
Every set of phase space points characterizes a face of W, with dimension one less than the 
number iV of points. So the types of faces of the kernel polytope, surfacing in a face of W, 
correspond to equivalence classes of sets of phase space points. For N = 7 there are five 
vertices pu ne , given by the different lines formed by subsets of the 7 prescribed points. It is 
not difficult to classify: For N = 6 there are two types, one type giving faces containing three 
vertices pn ne , the other two. For N = 5 there are two types, one containing two vertices, 
the other only one. For N = 4 either one vertex or none is present; iV = 3 is either a phase 
space line, giving one vertex at the symmetry center, or another triple with no kernel vertex. 
Also the edges, N = 2, contain no vertex. 

Note that the facts about isotropic states of qutrits, and the witness hyperplanes B Ptq which 
we use, as above in Theorem HJ is found by us in a new way as a byproduct of our special 
methods. See the next two Sect. IV Bl and Sect. IV CI where we proceed to find out more 
about the border between the separable and the entangled states. 



B. PT of our 9x9 matrices 

We represent the density matrices in the basis of product vectors 

\s-£,s) = \s-£)®\s) (39) 

and order them into groups of three, according to I. Inside each group we order according 
to s. So the global Hilbert space is represented as a direct sum 

Ci ©Cf =1 ©Ci 2 . (40) 

The projectors 

3 s,t 

do not mix the subspaces. A general matrix of W splits therefore into the direct sum 



P M = \ E u,fc( *~ t) I s -t,s)(t-£,t\ (41) 



P 



k,i \ k J \ k J \ k J 



with the 3x3 matrices 



,,111\ . / 1 w* w \ , / 1 WW* 

P = - 1 1 1 , P 1 = - I w 1 w* ), P 2 = - \w* 1 w I (43) 
° 1 1 1 1 / 6 \ w* w 1 / 6 \ w w* 1 
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Partial transformation is now the linear mapping 

\s-£,s)(t-i,t\ i-> \t-£,s)(s-£,t\ = \m-s,s)(m-t,t\ (44) 

with m = s +t — £ . 

It produces a new grouping of the basis vectors \m — s, s), according to m, and a new splitting 
of the global Hilbert space as 



"•3 rr\ rn3 rr\ ff~*3 
'm=0 



Ci =1 © Ci =2 . (45) 
The most general element of W is 

p = A e=0 © A e=1 © A i=2 (46) 

with 

j / d a* a \ 

A = - a d a* , tZ/ = ^ c M , a £ = w* CJM . (47) 
\a* a d I k k 

Partial transposition maps, as is demonstrated in the appendix, such a matrix p into 

, / d a 2 a* \ 

B®B®B, B = - a* d 1 a (48) 

\ai al d 2 J 

with three times the same 3x3 matrix. Thus a test for PPT of p reduces to a check for 
positivity of the matrix B. 

Now we apply the method and use the Peres criterion Ref. (l6j |: PPT, the positivity under 
partial transposition, is a necessary condition for separability. So NPT, non-positivity under 
partial transposition, implies entanglement. The missing detail in the proof of Theorem 
that the 72 triangular faces of W not corresponding to phase space lines contain no separable 
point, is contained in the following 

THEOREM 5 Consider a five dimensional face F of W, opposite to a triangular face 
which contains a separable point. F is spanned by the six Bell type states which are not 
located on the phase space line giving the separable state in the triangular face. The only 
entangled states in F, including its bordering faces, are two pu ne , and the edge joining them. 

PROOF By equivalence, we may assume that it is the vertical line with £ = 2, which gives 
the separable state in the triangular face, and which stays empty in F = {p = J2k c fc,o-ffc,o + 
Ylk c k,iPk,i}- The face F contains pu ne i= and pu ne i=x, and the edge joining them {ap Hnel=0 + 
(1 - a)pi inei=1 }. 

The general matrix, see ()46|) . in F is p = Ai=o@ At=\ ©0. It is transformed by PT, see 
to three times 

l ( d a\ 
B = o rf i a o 

This is only then positive, if ao = a\ = 0. That, in turn, implies, by (jUj), co,o = c^o = C2,o 
a, and cq,i = c^i = C2,i = 1 — ct. All other p are NPT, hence entangled. 
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The cases of centers of triangular faces not corresponding to a phase space line are 
represented by c 0i o = Ci j0 = c 0; i = |, with non-vanishing a and a^.D 

In geometric terms, this theorem is a statement about the N — 1 dimensional faces 
of W, spanned by N vertices with Bell type states. Up to iV = 4 there appear at most single 
separable points. For N = 5 and N = 6 we have to distinguish the equivalence classes. 
Already for N = 5, in the class not contained in a face F treated above, there appear new 
problems: Around the center, the state with coefficients co,o = C\ t o = C2,o = co,i = co,2 = §, 
there is a full four-dimensional set of states, which are outside the kernel polytope, but 
PPT. This can be observed by looking at the matrix B, obtained by PT of this center. It 
has the coefficients do = |, a = 0, a\ = a 2 = d\ = d 2 = |. These coefficients allow for 
small variations, without destructing the positivity of B. 

In the next application, looking into the interior of W, we study entanglement of mixtures 
of two orthogonal Bell type states and u. Note the generality of this case: We use the 
methods developed for W, but, as it follows from our Theorem Q we can choose any pair 
of mutually orthogonal Bell states, without reference to any special version of W. To apply 
our methods, we represent the state as 

p = 1 - ^ + ® 1 + aP lfi + PP 2>0 , with {l + 8a-8/3>0,l-a + 8/3>0}. (49) 
9 

This gives, besides a x = a 2 = 0, the no n- vanishing matrix elements 

l+2(a + /3) 



d 

d\ = G?2 



3 

1 - (a + (3) 



a = —+t—{(3-a). 



Now, p is PPT, iff B > 



& 0<di-\a \ 2 = 1 -(l-2(a + (3)-^(a + (3) 2 ^j - - a) 2 . (50) 

This describes, when the inequality is replaced by an equality, a nonlinear border between 
PPT and NPT states. Special points on this border are: 

• isotropic states, (3 = 0, border point at a = | 

• middle line, a = /3, border point at a + /3 = | 

Now, PPT is, for qutrits, no longer sufficient to guarantee separability, see Refs.[i3,E3,E3; 
"bound entanglement" may occur. So we have to use more specialized methods, to study 
the faces not covered by the Theorem and to analyze the interior of W. 

C. Constructions of witnesses 

An entanglement witness, EW P , gives a criterion, to show that a certain state with density 
matrix p is not contained in SEP, the set of separable states; see Ref. jiol ]. 

{EW P } = {K = | V<7 G SEP : Tr(aK) > 0, Tr(pK) < 0} . (51) 
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In this paper we are mostly interested in the structure of SEP. It is a convex set, and as 
such completely characterized by its tangential hyperplanes. The tangents itself are at the 
border of a larger set of hyperplanes which do not cut SEP. So we extend the meaning of 
witness and define SW, the set of structural witnesses: 

SW = {K = ^ | V<7 G SEP : Tr(aK) > 0, } . (52) 

Similarly, we define TW P , the set of tangential witnesses for a state on the surface of SEP 

TW P = {K = K ] ^ | V<7 G SEP : Tr{aK) > 0, Tr(pK) = 0} . (53) 

The set SW is convex and closed. It is also a linear cone: K G SW, a > 0, =>- aK G 
SW. Therefore the bounded set {K G SW, TrK'K < 1} contains all the mathematical 
information about SW; especially, that its boundary, that is TW, is closed. Also, that 
for every K the family {aK — (1 — a)u} has to intersect the boundary at some TW P . In 
geometric terms: In each family of parallel hyperplanes, in the Euclidean space of hemitean 
matrices, there are two tangential planes. Moreover, because of convexity, closedness and 
boundedness of SEP: For every K in the boundary of SW, there exists at least one p G SEP, 
such that K is a TW P . The boundary of SEP is our object of main interest. 
Here we analyze SEPfl W, and the symmetries of W are an important tool. That a symmetry 
of a state is reflected in symmetries of its witnesses seems intuitively clear, and has been 
used already, Ref. 0. We use this correspondence of symmetries in several details, so we 
formulate it explicitly: 

THEOREM 6 Consider a symmetry group Q , implemented by unitary and/ or antiunitary 
operators V g . Suppose that p is Q-invariant, t.i. \/g V g pV g x = p. 

1. If p is entangled, there exists a Q-invariant EW p . 

2. If p is on the surface of SEP, there exists a Q-invariant TW P . 

3. The subset of Q -invariant elements of SEP is completely characterized by the subset of 
Q-invariant elements on the surface of SW, which is the set of Q-invariant TW P . 

4- All the same is true, when SEP is replaced by the set of PPT- states, entanglement by 
NPT. 

PROOF We use the symmetrizing twirl operation, see Ref. |4|, K \— > (K)g. Here we use 
only finite discrete groups, so the Haar-measure is just summation, and 

9 

where \Q\ is the number of elements in Q. For every invariant p we have 

Tr(Kp)=Tr(K(p) g ) = Tr((K) g p). 

If K is an EW P , then also (K)g is an EW P , proving (1). If p is on the surface of SEP, there 
exists a TW P , say K, such that Tr(Kp) = 0. Then also (K)g is a TW P , proving (2). The 
set SW is convex and spanned by all convex combinations of TW p . The same is true for 
the subset of invariant structural witnesses; they are spanned by all convex combinations 
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of invariant tangential witnesses, on the surface. The invariant density matrices form also 
a closed convex set, of lower dimension. It is completely characterized by its tangential 
hyperplanes in this lower dimensional space. These are given through the invariant SWs. 
This proves (3). To prove (4), one observes that the PPT states form a closed convex set, 
and that the distinction between PPT and NPT is invariant under unitary and antiunitary 
transformations. □ 

For applications using the symmetries inside of W, we use later also the phase space 
reflections, implemented by local antiunitaries. So we had to consider also this kind of 
group action. 

As first application we consider the group of unitaries U^ i = 2P kj £ ~ 1 an d their products. 
They are reflections, U% e = 1. Our object W is pointwise invariant under this group, and 
its linear span is the largest set with this property. So, to study witnesses characterizing W, 
we have to consider the invariant opterator 

K = J2 «M P M • (55) 

This is sufficient to obtain all facts about the structure of SEP and PPT in W. Note: K is 

an EW P for some state, iff at least one K/~/ < 0. 

Now, once more, we use the "magic" of Bell type states. 

THEOREM 7 The operator 

K = J2 «M^M (56) 
is a structural witness i/f V0 G C 3 the operator 

M = X)«mWmI0><0|Wm (57) 

is not negative. 

K is moreover a TW for some p e W, iff 30, such that det = 0. 

PROOF Each separable state is a convex combination of pure product states \ip,r])(ip,r]\. 
So K is a SW, iff 

V-0,77 (-0,7717^1^,77) > 0. 
Now we use the definitions defined in Sect. ITD 

Pk,i=\^2w Kt \s,s){ft\W^ 

s,t 

and calculate 

(<ip,r}\P k/ \ip,r}) = 

s,t 
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giving 

U,, V \K\il>,Tl) = \(1>\MM- ( 58 ) 

Here we defined the vector G C 3 as |0) = ^2 s {v\ s )\ s ) , with the complex conjugated 
expansion coefficients of rj. What is true for all is then obviously true for all rj, and vice 
versa. 

If det = 0, there exists an eigenvector if) with eigenvalue 0. So K is TW P for a density 
matrix we can define explicitly by 

p = (\^,v)(^,v\)s- 

We use here the group Q with the Uk,e which we used in the first application of Theorem 
so p G W. On the other hand, if 3p such that K is TW P , one may expand p = (\if>, rj) (if), rj\ + 
Prest, with p rest also being separable. Then 

TiKp = -(iPlM^ip) + TrK Prest = 0. 

Each of the contributions has to vanish; so if) is an eigenvector of with eigenvalue 0, and 
det M = 0. □ 

First application of the Theorem [7| The well known optimal EW for a Bell type 
state. Consider P ,o- We use the symmetry of the phase space sub-group where we fix one 
point, e.g. (k,£) = (0,0), and mix all the other phase space points. An invariant witness 
has to have the form K = 7Po,o + (31. This gives 

M^ = 7 |0)(0| + 3/?||0|| 2 l. (59) 

We have used the representation of the unit operator on the global Hilbert space as 1 = 
J^ke-Pk/- This gi yes then as contribution to on C 3 the operator e Wk,e\4>) {(f>\W^ . 
This operator is invariant under the Weyl group, and its trace is 9||0|| 2 . This can only give 
3||0|| 2 1, as the contribution to M<f,. The eigenvalues of M^, for normed 0, are 7 + 3/3, (5,(5. 
So, if 7 = —3(5, K is the isotropic witness TW P for 

P = \Po,o + \u ■ (60) 

In this determination of EW, the choice of was completely irrelevant. This is connected 
with the high symmetry of Po,o- 

In a next application we consider fewer symmetries. We use the same methods as in Sect. IIVI 
in the equations from (JTTJ) to ()15j). Let Q, implemented by local unitaries or antiunitaries 
VgVg 1 , be the invariance group for p, an element of W. Choosing an invariant witness K, 
associated to the set of matrices M^, then every is unitarily equivalent to M x , with 
X = V~ x <f). This is seen by applying V g V~ l from left and its inverse from the right onto K 
in equation (|58jl . and calculating its action onto the matrix M^. 

For states and their witnesses which are located on a line in phase space, with an eventual 
part proportional to w or 1, this brings essential simplification. By equivalence, we may 
consider the line {(0,0) . . . (2,0)}. All the Wk,o, consequently the p^o, and of course also 
the unit operator, are invariant under the group of unitaries ^ s e tS ^\s)(s\. The consequence 
is, that each is equivalent to with real valued non negative vector \<p\. 
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For general witnesses there remain the phase space translations V g , combined with V g as 
local unitary operators. They act onto the Weyl operators by multiplication with phase 
factors which cancel in the action onto K. The consequence for equivalences of are 
the symmetries of detM^ under cyclic permutation (p s hh> <fi s+ i, and the discrete phase 
twirling <ft s i— > w s <p s . One knows therefore, that the determinant depends on the <j) s and their 
conjugates in the form of symmetric polynomials, invariant under the discrete phase twirling. 
This allows permutation-symmetric sums with contributions from \(j) s \ 2 , from 0q0*02, etc. 
But it forbids contributions as 4> 2 S , ^\^\<^\ etc. 

Combining these results, it is not difficult to calculate the determinants for witnesses located 
on the phase space line i = 0, mixed with 1: Consider 

K = 4l + ^7fcP fc ,o> ( 61 ) 

k 

related, when \\(j>\\ = 1 to the matrices = Al + ^ k jkWkfl\(j))((j)\W^Q. The matrix is 
written explicitly in the appendix. The determinant works out as 

detM^A 3 + ||0|| 2 (7o + 7i + 7 2 )A 2 (62) 
+ 3 (|0 o | 2 |0i| 2 + |0i| 2 |02| 2 + |02| 2 |0o| 2 ) (7o7i + 7l72 + 727o) A 
+ 27|0 o | 2 |0i| 2 |0 2 | 2 7o7i72. 

The analysis follows in the next sections. 



D. Some details in the structure, analytical 

The strategy for the exploration of the structure of SEP is to find the set of tangential 
witnesses as follows: Analyse the operators K = t Kk,ePk,e by way of studying the set of 
matrices associated to each single K. If these matrices are positive for all 0, then K is a 
SW . If there is a 0, such that det = 0, then K is a TW. If one has "enough" TW p s, one 
can determine the {/?} in the boundary of SEP. Consequently, they obey TrKp = 0. How 
many of these witnesses are "enough" , depends on the symmetry of the subset of states one 
is studying. High symmetry restricts and simplifies the study. 

We study the subset of states with components located on a phase space line, mixed with 
uj. By equivalence, it is sufficient to study one special line. We choose that with t = 0. The 
states are p = a J2k K fc,o-Pfc,o + (1 — ct)uj with J2k K k,o = 1- Each of these states is invariant 
under vertical shear and horizontal reflection. This symmetry group implies that we may 
restrict the study of witnesses to 

K = \^l + J2^Pk,o, (63) 

k 

as at the end of Sect. IV CI Note that the parameter A cannot be negative for witnesses, 
since Tr Kpu ne ^ = A/3, but Tri^cr should be non-negative for separable states. And we 
know already, that all the pu ne are separable, Theorem 

Especially simple are those operators, where A = 0: K is a TW, iff all > 0. Because, if 
all the factors 7^ are non-zero, then each is a sum of positive operators. But, if one of 
them is negative, then, with = J2 S 75 I s )' ma t rrx has a negative eigenvalue. Such 
a K with non-negative 7^ is no EW, but tangential to the face of W spanned by the Pk,i^o- 
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FIG. 2: (Color online.) Here the space of the density matrices p = ~q w + ai-bo + f (-P10+-F20) is 
shown. The green triangle shows the positivity condition. The blue lines are the inner fence (kernel 
polytope) and the outer fence (enclosure polytope) for the boundary of separability. The dotted 
area shows the values of a and (3 which are positive under PT. Here PPT equals separability. 



In the study of the operators with A > 0, it is enough to consider A = 1, because the 
witnesses form a cone, all parameters may be scaled. The investigation, whether > 0, 
is now done by investigation of the characteristic polynomial S c (/i) = det(M^ — ii). Since 
Mtf, is hermitian, this polynomial, S c (/i) = n^/ 1 ?' — a0> nas om y rea ^ zeroes /ij. We have to 
demand that they are not negative, and this is the case if and only if 

• the second derivative of S c at fi = is not negative, 

• The first derivative there is not positive, 

• H c (yU = 0) is not negative. 

These are conditions for SW. To get TW, we need one eigenvalue equal to zero, and 
strengthen the last condition to 

• H c (// = 0) = 0. 

The characteristic polynomial is given by replacing A with 1 — /x in the equation (J62j) . We 
use real valued and the abbreviations 

A := 7o7i + 7i72 + 727o, B := 707172 (RA \ 

f A := 3(06^ + 0? 2 + 2 2 }) f B ;= 270 2 2 2 • 
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The conditions for a tangential witness are given by calculating the derivative S c , using 
11011 = 1, 

• 3 + 2^ fe7fc + min (A-/ A ) > 0, 

• 1 + E fc 7fc + min (A • f A + 5 • / fl ) = 0. 

The minima over normalized wavefunctions are evaluated in the appendix: 

mm(A-f A ) = min{0,v4}, (65) 

min(A .f A + B-f B ) = min{0, -A, A + B}. (66) 

<p 4 

Because of f!65|) . the first of the conditions as stated above for witnesses is weaker then 
the second one. And the parameters for TWs have to fulfill only one inequality and one 
equation: 

3 + 2 7o + 4 7 + min{0, A} > (67) 
1 + 7o + 2 7 + min{0, -A, A + B} = ; (68) 

where we use now the parameters 

7 = ^(7i + 72), $ = -(71-72). ( 69 ) 

Using them we get 

A = 27q7 + 7 2 - S 2 , A + B = 2 7o 7 + 7 2 + 707 2 . (70) 

In the search for TW P for p symmetric under vertical reflection, i.e. K10 = ^2,0 5 one can 
restrict the search to operators K which have the same symmetry, that is, they have 5 = 0. 
We explore the set of witnesses starting from the isotropic witness, with A = 1, 70 = — 1, 
7 = 5 = 0. First we list all those TWs one gets, then we indicate the "proof". 
In the set of results we find four distinguished regions for the parameters: 



a) A = 


1. 


7>0, 


7o = 


-i; 


b) A = 


1. 


0> 7 >-§, 


7o = 


-1 - 27 < 


c) A = 


1, 


1 = 'I 


7o > 


1 . 

3 ' 


d) A = 


0. 


7>0, 


7o = 


1 - 7 > 0. 



In the parameter region a) one has A = / ~f 2 — 27 > —1, so the l.h.s. of Eq. (JfiTj) is positive; 
and min{0, \A, A + B} = —27, so Eq. ()68)1 is true. At one end of the region, i.e. in the limit 
7 — > 00, one may rescale the parameters and observe that they approach A = 7 = 0,7>0, 
one end of region d). At the other end of region a), which is also the beginning of region 
b), the common witness at the edge of these regions is the isotropic witness with 7 = 
and 70 = —1, corresponding to the hyperplane -E>o,o, defined in Eq. flBlIj) and in Eq. (JMj) . In 
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the parameter region b), succeeding a), neither A nor A + B is negative, since 7 < 0, as 

long as 7 > — |, with 7 related to 7 by (jSBj)- in the succeeding region c) one has A < 0, 

and ^A < A + B. At the end of this region rescaling leads here, in the limit 70 — > 00, to 

A = 7 = 0, 70 > 0, this is the other end of region d). The round trip is finished. 

For regions a), b) and c), we get EWs, except for 7 = 70 = — | in region b). There we get 

the TW p for all the p in the triangular face of W with the P^.o a t the vertices. 

The linearity of the relations implies, that in each regional set of witnesses each A' is a TW P 

for (at least) one common p. For each K at and endpoint of a region, there exists a linear 

face of SEP, for which K is tangential j2l|. 

The vertex points of SEP, corresponding to the linear regions of the witness-parameters, 
are: 

a) O a = L) + §-Po,Q ~ gPlfi ~~ 5-^2,0) 

b) Ob = Plinei=0] 

c) a c =l(u- §P ,o + |Pi,o + |P 2 ,o); 

d) CT d = ^(u; — Pline£=o) = \{Plinel=l + P\,inei=2)- 

The hyperplane -Bo,o> corresponding to the isotropic witness, contains four pu ne , associated 
with the four phase-space-lines through the phase-space point P 00 . One of them is <r b ; the 
other three have a a in their middle: 

Va = \ Q(P ,Q + P ,l + P , 2 ) + ^(P ,0 + Pl,l + P 2)2 ) + ^(P ,0 + ^2,1 + P, 2 ) 

The part of the boundary of SEP between cjb and a c is outside the kernel polytope and inside 
the enclosure polytope. It crosses the ray from uj to |(Pi,o+p2,o) at v-mid = l^cu+P^o+^o)- 

E. Some details in the structure, numerical 

Our method of numerical analysis is a variation of the strategy we use in the previous section. 
There we calculate first the whole set of tangential witnesses K, then we find the p on the 
border of SEP via the condition Tr(pK) = 0. Here we use this condition from the beginning 
and reduce in this way the set of parameters which have to be varied. We find that explicitly 
in the following way. 

Let us here consider again the density matrix p = ' UJ + a p l +[3p 2 where the state space 
is visualized in Fig. El Via Theorem[7|the operator K = t + aPoo + ^-Po + CP20 is a structural 
witness iff the matrix = 1 + aWo t0 \4>) (4>\W 0)0 + &Wi,o|0)(0|Wi 1 o + cW 2 ,o|0)(0|W2,o is non 
negative V0. We are moreover interested in the tangential witness, therefore we search for 
Tr(pK) = which leads to a = 3 + b ( 1 + 8a ~^+^( 1 ~ Q! + 8 / 3 ) _ Consequently, we have to search for 
the border where for a given a and the variation over b, c and one eigenvalue of changes 
from negative to positive (the two others are positive). The found minimal (3 characterizes 
then the border state for which a separable state changes to an entangled state, and K is the 
tangential entanglement witness. We did not include any further symmetry constraints into 
the calculation and found that the analytical symmetry results as described in the previous 
section and in the appendix are confirmed, e.g. we have only to vary over real vectors <fi. 
In the region where no bound entanglement was found, the optimized parameter (3 agreed 
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FIG. 3: (Color online.) Here the space of the density matrices p = ~q uj + aP\o + PP20 is 
shown. The green triangle shows the positivity condition. The axes joining uj and P10 or P20 
represent the isotropic states. The blue lines are the inner fence (kernel polytope) and the outer 
fence (enclosure polytope) for the boundary of separability. The dotted area shows the region of a 
and (3 where p is PPT. For both parameters a, [3 > PPT gives the boundary between separable 
and entangled states. For a or (3 negative we find bound entanglement which can be seen in the 
enlarged picture. The dashed curve shows the border of PPT and for the points the witnesses were 
explicitly numerically calculated, see Sect. IV El 

with the PPT calculated f3 numerically up to 10~ 8 . The largest difference between the PPT 
boundary and the separability boundary is of the order of 10~ 2 and decreases to zero for a 
or (3 approaching 0, the isotropic state, see Fig. El 

VI. SUMMARY AND CONCLUSIONS 

We consider the state space of two qutrits where we restrict ourselves to locally maximal 
mixed states. Whereas for qubits every locally maximally mixed state can be diagonalized 
by the magic Bell states and therefore embedded into a magical tetrahedron, this is not 
true for qutrits. However, we show that a kind of analogue to the magic tetrahedron can be 
denned for qutrits: the magical simplex W. 

Starting from a certain maximally entangled pure state, a Bell type state, we obtain by 
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applying only on Alice side the Weyl operators nine other orthogonal Bell type states. The 
Weyl operators are used to describe the discrete classical phase space. This discrete classical 
phase space representing the algebraic relations of the Weyl operators enables us to describe 
the local transformations of the state space of interest and are very useful for several proofs 
in this paper. The mixtures of all Bell type states form the simplex W which is then the 
main object of our investigations. 

We show explicitly how to construct a version of W. A certain version is fixed by defining 
3 Bell type states. The simplex W can be embedded in a 9-dimensional Euclidean space 
equipped with a Hilbert-Schmidt norm and an inner product. 

Transformations of W onto itself can be considered as transformations of the discrete classical 
phase space. Thus the symmetries and equivalences can be studied by this means. 
Then we investigate the question of the geometry of separability. We start with construct- 
ing two polytopes, an inner (kernel polytope) and an outer (enclosure polytope) fence for 
separability. They define linear entanglement witnesses but are in general not optimal. The 
outer fence, the closure polytope, has the same geometric symmetry as W. 
Hereupon we explicitly study two representative cases. We consider the state space of all 
density matrices which are mixtures of the total mixed state and two Bell states. We apply 
the partial positive transformation on one subsystem (PPT) which detects entanglement. 
The obtained witness is no longer a linear one. We show how entanglement witnesses 
can be constructed and apply it to the density matrices under consideration. We find after 
optimizing the entanglement witness by analytical and independently by numerical methods 
that there is indeed bound entanglement for negative mixtures of one of the Bell states. The 
result is also visualized in Fig. El 

The second case we study is the state space of density matrices which are mixtures of the 
total mixed state, one Bell state and an equal mixture of two other Bell states. We find that 
it has only linear entanglement witnesses and that no bound entanglement can be found, 
visualized in Fig. 

Summarizing, we could give a full geometric structure of the subset of bipartite qutrits under 
investigation. We think that this will help to find a good characterization of the whole state 
space and to investigate measures for entanglement for higher dimensional systems. 
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VII. APPENDIX 

1. A physical model 

Each party has a system consisting of a ring shaped molecule. In this ring there are three 
symmetric located possible places for a single itinerant particle. Locating this particle at any 
of these places corresponds to our three basis vectors \s). In the entangled state, described 
by the vector f^, the index I denotes the angular correlations between Alice's and Bob's 
particles. Concerning measuring of locations, for I = they are to be measured at the places 
at the same angles. For the two other cases, Alice's particle is rotated relative to Bob's. 
So the number £ is the quantum number for the observable SBob — s Alice] an d the index 
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k is the quantum number for the total angular momentum (component orthogonal to the 
rings). Again "—1"= 2, due to the finiteness of the system. These two operators commute, 
while their individual contributions from one party do not; comparable to the commuting 
of xb — %a with pb + Pa- The set of their eigenvalue pairs (k,£) is the discrete classical 
phase space. 

2. Examples of maximally mixed states which do not fit into W 

An example for the observation (JTJ) in Sect. [TTT| Define p = ||^ r )(^ r | + §|$)($|, with = 
|0) (g> |0), and |$) = ^(|1) <8> |1) + |2) <8> |2)). This density matrix is in a unique way 
diagonalized, - with non Bell states. 

As an example for the observation (J2J) we consider p = ^2 c a \^f a ) (& a \, with three different 
orthonomal normalized Bell vectors ty Q and three different expansion factors c a . Since such 
an expansion is just the expansion into projectors onto eigenvectors, it is unique. So, if p 
can be embedded into W (or a unitary equivalent), the expansion must be an expansion 
into the P^ (or into a unitarily equivalent set). 

Now we give an example of three Bell type projectors which cannot together be embedded 
into W: Take two of the projectors as Po,o and Pi,o, the third one as |$)($|, with |$) = 
75 J2 S w2s ^(^\ s ) + 2(s — 1) — \s + 1)) ® \s). This is a Bell vector, orthogonal to fio.o and f^o- 
But it is not orthogonal to f^O; so I^K^I cannot be any P^. Now try a transformation and 
embed the first two projectors as Po,o — P'ht anc ^ ^i,o = P'k+mi+n a unitary equivalent 
version of W. Consider the mapping between these projectors by Weyl operators. As the 
following equation shows, they are fixed up to a phase 

Pl,0 = PL+mJ+n = (U® l)P' k/ {U ® l)t = {U® l)P Ol0 (C/ t ® 1) ■ (71) 

Taking the matrix elements with \s) <8> \t), identifying U with W mn and comparing with the 
relation between Pq^ and P^o gives the equations 

(s\W^ n \t) = (s\U\t) = e^(s\W lfl \t) = eV6 s , t . (72) 

And this implies, by applying W mn once more, that P'k+2mi+2n = -^2,0 > an d the third Bell 
state does not fit into the equivalent version of W either. 

An example for the observation (JHJ): Take the three Bell states Qk,2 out of W and replace 
them by 

®k = ^^2 wkSa s\s-€) ® \s). 

s 

They are orthogonal, span the same subspace as the deleted flk,2, but define other states, 
unless the phase factors a s are chosen in a very special way. Together with the remaining 
Qk,o and f^i they form a complete set of orthogonal Bell vectors, but nothing equivalent to 

w. 

3. Matrices representing the Weyl operators 

The basis vectors are 
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The Weyl operators W k ^ arranged according to the appearance of the indices in the phase 
space are 



k = 2 



k = 1 



k = 



1 

w* 
w 

/ 1 
w 
\0 






(74) 



Complex conjugation interchanges the lines k — 2 and k = 1. 

The transformation producing unitary operators are represented as 




1 

U v = ( 1 | . (75) 
u>* 



^. TTte partial transposition 

The basis elements are the product vectors \s — £,s). They are arranged in groups of three, 
according to £. Inside each group the ordering is according to s. The index-pairs (s — £, s) 
denoting the rows are written on the left side. Partial transposition induces a new splitting 
of the Hilbert space in subspaces. They are characterized by m, when we write the index- 
pairs now as (m — s, s). We mark the different m by different typefaces; a for m = 0, a for 
m — 1, a for m — 2. (But the numbers are the same, independent of the typeface!) 
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(77) 



5. The matrix for witnesses on a line, mixed with the unit 

We use real valued S . This is possible because of the invariances as described at the end of 
the Sect. EH] 

A + 0o(7o + 7i + 72) 0o0i(7o + ^*7i + W72) o 02(7o + W7i + w *l2) \ 
0i0o(7o + wr/i + w*72) A + 0^(7o + 7i +72) 0i02(7o + ^*7i + wry 2 ) . (78) 
0200(70 + W*7i + W7 2 ) 0201 (7o + W7i + w *l2) A + 0^(70 + 71 +7 2 ) J 

6. The minima for the functions of <f>, used in Sect. I V LA 
We use real valued, normalized wavefunctions and describe them with the two parameters 

z = <& G[0,1], (79) 

* = \{ti-ti) (80) 

The minima of the functions of can be evaluated as minima of functions of z and x in 
the triangular region determined in (fTD~J) and (fHUJ) . The functions defined in the Sect. IV Dl 
in Eq.([64p are 

f a = ^(l + 2z-3z 2 )-3x 2 , (81) 
f B = 27 'z - z 2 ) . (82) 

Because of the terms quadratic in x, the extrema are attained either at the boundary of the 
triangle, or at the line x — 0. One finds the minimum of f A at the triangle-vertices, equal 
to 0. The maximum is attained in the center, at z = |, equal to 1. Now one has to respect 
the sign of A: 

mm(A ■ f A ) = min{0, A} . (83) 



In the combinations of f A with also the maximum of at the border of the triangle 
comes into consideration. It is found at the center of a border line and it is equal to |. 
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The non negative function /# is zero along the whole boundary and has its maximum, equal 
to 1, at the center. The minimum of Af A + -B/b is zero, attained at a vertex, if both A and 
B are non-negative. It is attained at the center, and equal to A + B, if both A and B are 
non-positive. For different signs of A and B one has to analyse Af A + Bfs along the line 
x = 0. Consider 

/a + Cf B = - (1 + 3z + 9Cz(l - z)) 

as functions of z, depending on the parameter C = B/A. For every C they have fixed values 
at z — - there they are | - and at z = 1 - where they are 0. For every C the first 
derivative at z — | is zero. There is either a local minimum - if C < — | -, a local maximum 
- if C > — | -, or a saddle point. Now we consider the function with special values of the 
parameter C: For C = —1, it is zero at z — 1/3, which is a local minimum. For C = —1/4, 
the local maximum at z — 1/3 is equal to the maximal value | at the border, at z — 0. 
Inside the region < z < 1 the set of functions is pointwise monotone increasing in the 
parameter C. So for C in between the special values, the minimum is zero, the maximum is 
|, both attained at the border. For C < — 1, the minimum, 1 + C, is attained at the center, 
the maximum, |, at the border. For C > — |, the minimum, zero, is attained at a vertex, 
the maximum, 1 + C, at the center. 

The result can be stated in a unified way as the equation 

min(A -f A + B-f B ) = min{0, jU, A + 5} . (84) 

x,z 4 

It means, that the search for the minimum over all (f) can be reduced to considering only 
three different vectors: 

£(!)■ Mi)- (85) 
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